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Abstract 

The  boundary  value  problem  of  sound  generated  at  an 
oscillating  wall  propagating  into  a  half  space  is  solved 
for  a  linear  Boltzmann  equation  with  a  general  cutoff 
Intermolecular  potential  and  an  arbitrary  boundary  condi- 
tion at  the  wall.   For  a  small  bounded  domain  (two  walls 
closer  than  a  mean  free  path)  proof  of  the  existence  of  a 
solution  is  relatively  simple.   In  a  half  space  the  solu- 
tion is  proved  to  exist  for  sufficiently  high  frequency 

00  >  o)   (o)  is  comparable  to  the  mean  collision  frequency). 
^   o  ^  o 

Asymptotic  expressions  show  the  disappearance  of  a  conven- 
tional sound  wave  at  high  frequency.   This  is  consistent 
with  known  results  for  a  relaxation  model  of  the  Boltzmann 
equation  and  qualitative  estimates  for  the  actual  Boltzmann 
equation  based  on  the  dominance  of  the  continuous  spectrum 
at  high  frequency.   At  a  fixed  arbitrary  distance  from  the 
wall,  the  dominant  high  frequency  disturbance  consists  of 
just  those  particles  emitted  by  the  wall,  decreased  by 
scattering,  but  without  any  compensating  contribution  from 
the  creation  term  of  the  collision  operator. 
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1 .   Introduction 

The  propagation  of  sound  waves  from  an  oscillating 
wall  Into  a  gas  has  been  treated  by  many  approximate 
methods,  but  no  satisfactory  mathematical  solution  has 
been  given  for  the  Boltzmann  equation.   An  elegant  solution 
to  this  problem  for  a  model  of  the  Boltzmann  equation  has 
been  obtained  by  Weltzner  [1].   There  are  several  high 
order  moment  and  polynomial  approximations  [2],  [5],  [^]; 
but  it  is  difficult  to  assess  their  accuracy  since  they  are 
essentially  ad  hoc  and  insensitive  to  the  many  singular 
features  of  the  solution.   An  elaborate  model  of  the 
Boltzmann  equation  has  been  used  by  Slrovich  and  Thurber  [5] 
but  only  to  evaluate  normal  modes,  not  to  solve  the  boundary 
value  problem.   A  polynomial  approximation  in  the  neighbor- 
hood of  free  flow  has  been  described  which  is  much  more 
likely  to  give  an  accurate  representation  in  the  high 


1.  Harold  Weltzner,  "Steady-State  Oscillations  in  a  Gas," 
to  appear  in  Proceedings  of  the  Fourth  International 
Symposium  on  Rarefied  Gas  Dynamics,  Toronto,  1964. 

2.  C.  S.  Wang  Chang  and  G.  Uhlenbeck,  "Propagation  of  Sound 
in  Monatomic  Gases,"  Engineering  Research  Institute, 
University  of  Michigan,  1952. 

3.  L.  Slrovich,  Physics  of  Fluids  6,  2l8  (I965). 

4.  G.  L.  Pekeris,  Z.  Alterman,  L.  Flnkelsteln,  and  K. 
Frankowskl,  Physics  of  Fluids  5,  I608  (1962). 

5.  L.  Slrovich,  J.  K.  Thurber,  Journal  of  the  Acoustical 
Society  37,  329  (1965). 
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frequency  regime  [6];  but  again  the  error  Is  hard  to 
estimate  and  not  all  singular  features  are  visible. 

There  are  two  major  questions:   the  existence  of  a 
solution  to  the  problem,  a  decidedly  nontrlvlal  matter; 
and  descriptive  properties  of  the  solutions  which  are 
emphatically  singular,  especially  at  high  frequency. 

With  regard  to  existence,  there  has  been  a  substan- 
tial development  for  the  Initial  value  problem  (with 
trivial  boundaries)  [7],  even  to  some  extent  In  nonlinear 
problems  [8].   The  basic  tool  In  this  theory  Is  the  recog- 
nition that  a  certain  class  of  cutoff  Intermolecular 
potentials  yields  a  mathematically  accessible  Boltzmann 
equation.   But  certain  a  priori  estimates  which  are 
crucial  for  solution  In  the  large  In  the  Initial  value 
theory  are  not  available  In  boundary  value  problems.   For 
this  reason  only  existence  In  the  small,  e.g.,  between 
two  walls  which  are  sufficiently  close.  Is  an  easy  matter. 
We  use  the  high  frequency,  cjd  >  co  ,  to  Improve  our  estimates 
enough  to  extend  the  solution  to  an  entire  half  space. 


6.  D.  Kahn  and  D.  Mlntzner,  Physics  of  Fluids  8,  1090  (1965)' 

7.  Harold  Grad,  "Asymptotic  Theory  of  the  Boltzmann 
Equation,  II"  In  Rarefied  Gas  Dynamics,  ed.  J.  A. 
Laurmann,  Academic  Press  (1965)- 

8.  Harold  Grad,  "Asymptotic  Equivalence  of  the  Navler- 
Stokes  and  Nonlinear  Boltzmann  Equations,"  In  Proceedings 
of  Symposia  In  Applied  Mathematics,  Vol.  17,  American 
Mathematical  Society  (1965)- 
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But  the  question  of  solution  of  even  the  steady  flow  linear 
Boltzmann  equation  In  sizable  domains  is  largely  untouched. 
The  problem  of  making  plausible  rather  than  rigorous  asymp- 
totic estimates  (e.g.,  near  the  wall,  far  away,  at  very 
high  frequency,  for  low  or  high  speed  particles)  is  not 
so  difficult  and  is  not  significantly  harder  for  the  actual 
Boltzmann  equation  than  it  is  for  a  model  equation. 

The  main  qualitative  feature  of  all  problems  in  kinetic 
theory  is  the  transition  from  free  flow  to  continuum  flow. 
For  the  initial  value  problem  an  essentially  complete  answer 
is  known  [9],  [10].   For  a  very  short  time  there  is  an 
approximate  free  flow;  this  is  followed  by  a  readjustment 
to  an  approximately  local  Maxwelllan  on  a  time  scale 
comparable  to  the  mean  collision  time;  then  there  is  a  much 
longer  era  of  macroscopic  fluid  motion  and  an  ultimate 
decay  to  equilibrium.   Similarly,  in  a  steady  boundary 
value  problem  we  can  expect  something  like  free  flow 
(characterized  by  two  distinct  half -distributions)  very 
close  to  the  wall;  a  kinetic  boundary  layer  on  the  order 
of  the  mean  free  path;  and  a  subsequent  fluid  regime 
governed  by  macroscopic  equations.   But  instead  of  a  distant 
behavior  which  is  fluid-like,  there  is  (at  least  for 
moderate  frequencies)  only  an  intermediate  zone  in  which 


9.   Harold  Grad,  Physics  of  Fluids  6,  14?  (I963) 
10.   Harold  Grad,  Journal  SIAM  13,  259  (1965). 
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a  sound  wave  can  be  perceived.  This  Is  followed  by  a 
complex  kinetic  ultimate  decay  (discovered  for  the  model 
equation  by  Weltzner  [l]).  This  ultimate  kinetic  regime 
also  exists  In  the  Initial  value  problem  [7]  for  a  class 
of  soft  potentials  which  Is  unrealistic  in  neutral 
gases 5  in  the  boundary  value  problem  this  effect  Is  more 
c  ommon . 

The  reason  for  this  phenomenon  is  the  presence  of  a 
continuous  spectrum  extending  to  the  origin.   This  fact 
alone  provides  a  simple  qualitative  description  even  for 
the  actual  Boltzmann  equation  [11].   The  non-exponentlally 
decaying  contribution  from  the  continuum  ultimately 
dominates,  even  though  the  more  highly  excited  normal 
mode  is  initially  more  Important.   But  for  sufficiently 
high  frequency  the  normal  mode  (i.e.,  sound  wave)  is 
completely  submerged  at  any  distance  from  the  wall. 

In  a  boundary  value  problem  in  which  there  is  alge- 
braic rather  than  exponential  decay  at  infinity,  such  as 
flow  around  a  body,  the  effect  of  the  continuum  may  not 
be  evident  [ 12] . 

Another  source  of  unusual  and  singular  behavior  arises 


11.  Harold  Grad,  "Theory  of  the  Boltzmann  Equation,"  to 
appear  In  Proceedings  of  the  Eleventh  International 
Congress  of  Applied  Mechanics,  Munich,  19^4. 

12.  Harold  Grad,  "Flow  Equations  in  a  Rarefied  Atmosphere," 
in  Proceedings  of  Symposium  on  Aerodynamics  of  the 
Upper  Atmosphere,  Project  Rand  (1959)- 


-  4  - 


from  the  very  slow  and  very  fast  particles.   The  former 
are  disruptive  because  the  free  path  Is  arbitrarily 
small.   There  can  be  very  rapid  spatial  gradients,  and 
at  the  wall  itself  infinite  gradients  (see  Sec.  3).   This 
Is  perhaps  the  most  serious  hindrance  to  the  development 
of  a  general  existence  theory.   The  fast  particles  are 
characterized  by  unbounded  free  paths  (except  in  the 
special  case  of  hard  spheres),  and  they  dominate  the 
behavior  far  from  the  wall;  this  gives  a  more  intuitive 
picture  of  the  origin  of  the  continuous  spectrum.  An 
Important  quantitative  point,  that  the  collision  frequency 
is  velocity  dependent,  makes  it  quite  important  to  con- 
sider the  actual  Boltzmann  equation  rather  than  a  model. 


2.   The  Spectrum 

The  distinction  between  a  discrete  and  a  continuous 
spectrum  is  basically  mathematical,  but  the  qualitative 
difference  in  the  behavior  of  solutions  is  so  great  as 
to  yield  very  important  qualitative  physical  differences. 
In  a  normal  mode  all  the  essential  properties  of  a  wave 
belong  to  the  gas  alone  (e.g.,  the  phase  speed  and 
decrement);  the  initial  and  boundary  values  only  serve 
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to  determine  the  level  of  excitation.   But  a  mode  which 
results  from  a  continuum  never  loses  sight  of  the  entire 
initial  or  boundary  data;  It  does  not  have  a  definite 
speed  or  a  definite  wavelength,  nor  does  it  decay  exponen- 
tially.  And  each  of  these  properties  is  inextricably 
interwoven  in  characteristics  of  both  the  medium  and  the 
boundary.   Thus,  although  one  can  give  mathematical 
formulations  which  combine  both  possibilities  within  a 
common  framework,  it  is  necessary  to  distinguish  the  two 
in  order  to  understand  the  physical  consequences. 

There  are  two  sources  of  continuous  spectra  in  kinetic 
theory.   The  first  and  more  obvious  results  from  streaming. 
The  operator  ^  •  ^f/^x  has  the  entire  imaginary  axis  as 
its  spectrum.   In  the  form 

||.e.||=0  (2.1) 

(free  flow)  we  recognize  that  there  are  no  normal  modes. 
Yet  all  macroscopic  moments  decay  to  equilibrium  as 
t  — >oo,  not  at  a  predetermined  rate,  but  at  a  rate  which 
depends  on  the  initial  data  [13]. 

The  second  source  of  a  continuum  is  in  the  collision 
operator  L.   This  property  has  only  recently  been  recog- 
nized and  is  found  to  accompany  almost  all  intermolecular 


13.   Harold  Grad,  Communications  on  Pure  and  Applied 
Mathematics  Ik,    323  (1961). 
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potentials  [7].   In  the  form 

II  =  L[f]  (2.2) 

(spatially  -uniform)  we  observe  no  simple  exponential 
decay  to  equilibrium;  b\it  for  the  most  Interesting  case 
of  hard  potentials,  all  solutions  do  decay  faster  than  a 
given  exponential. 

The  spectrum  In  a  more  complicated  problem  reflects 
the  interaction  between  the  streaming  operator  and  the 
collision  operator.   We  write  the  full  linear  Boltzmann 
equation  in  the  form 

||.e-||  =  L[f]  (2.3) 

where 

L[f]  =  K[f]  -  vf  .  (2.4) 

L  is  the  linearized  collision  operator  and  is  non- 
positive;  the  origin  is  a  multiple  point  eigenvalue  In 
virtue  of  the  conservation  of  mass,  momentum,  and  energy. 

The  singular  part  of  L  is  given  in  terms  of  the  collision 

2     2     2     2 
frequency  v(^),  ^  =  ^  +  ^^  +  ^-   Po^  ^   certain  class 

of  repulsive  hard  force  laws  (this  includes  all  power 

laws  with  exponent  higher  than  the  Maxwelllan  fifth  power 

and  includes  hard  spheres),  v  is  monotone  increasing  in  ^ 

and  is  bounded  away  from  zero,  v(0)  >  0.   We  normalize 
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the  time  scale  such  that 

v(0)  =  1  ,  (2.5) 

thus  v(^)  >  1.   The  nonslngular  part,  K,  of  the  collision 
operator  is  a  bounded,  compact  integral  operator  [7]. 

From  these  properties  we  can  predict  the  location 
of  the  continuous  spectrum  (also  the  conservation  eigen- 
values), in  some  cases  as  a  consequence  of  proved 
theorems,  in  others  as  plausible  conjectures.   In  addi- 
tion to  free  streaming  (2.1)  and  spatial  uniformity 
(2.2),  we  consider  the  initial  value  problem  for  a  fixed 
wave  number  k, 

II  +  i(e-k)f  =  L[f]  ,  (2.6) 


steady  flow 


e  -11=  L[f]  ,  (2.7) 

and  separated  time 

icof  +  ^  -11  =  L[f]  .  (2.8) 

We  have  already  identified  the  continuum  in  free 
flow  as  the  entire  imaginary  axis  (real  cu) ,  Fig.  1. 
Under  spatial  uniformity  (2.2),  the  continuum  is  real 
and  consists  of  all  values  taken  by  -v(|).  Fig.  2.   The 
Initial  value  problem  (2.6)  has  as  its  continuum  the 
values  taken  by  -v(^)  -i^-k.  Fig.  3a  for  small  k  and 
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Pig.  3b  for  large  k.   The  perturbed  multiple  eigenvalues 


Fig.  1 


Fig.  2 


in  Fig.  3b  have  probably  disappeared  to  the  left  of  the 
continuum. 


Pig .  3a 


Fig.  3b 


In  a  steady  flow  (2.7),  the  continuum  covers  the 
entire  real  axis.  Fig.  4a,  since  It  takes  all  values 
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■5^ 


Fig.  4a 


Pig.  4b 


-v(^)/^-,  (t-.    Is  the  component  in  the  differentiated 
direction,  say  x.,  ) .   We  have  drawn  half  of  the  continuum 
dashed.   In  a  half  space,  boundedness  at  infinity  should 
eliminate  half  of  the  spectrum.   The  point  eigenvalues  at 
the  origin  are  not  perturbed  and  correspond  to  three  exact 
solutions:   f  =  const.,  T  =  linear,  v  =  linear,  corresponding 
to  a  trivial  solution,  simple  heat  flow,  and  simple  shear 
flow  (the  first  Chapman -Enskog  or  Navier-Stokes  term  is 
exact).   For  the  case  of  hard  spheres  only,  v{^)/^      is 
bounded  away  from  the  origin,  and  the  spectrum  is  as  shown 
in  Pig.  4b. 

With  time  separated  (2.8),  the  steady  flow  spectrum 
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becomes  complex.  Pig.  5a  for  small  co  and  Fig.  5t>  for 
large  cd;  for  hard  spheres  the  spectrum  Is  as  In  Fig.  3 
(plus  a  mirror  Image  spectrum  In  the  right  half  plane). 


o 


\ 


Pig.  5a 


Pig.  5b 
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2     2     2 
Since  v  is  a  function  of  2,.,  +  f.^  -^  ^^  and  only  the 

component  i-.    multiplies  Bf/^x-,  ,  the  complex  spectra  are 

not  curves  but  entire  regions.   The  spectrum  covers  the 

shaded  areas  in  Pigs.  3  and  5- 

For  a  relaxation  model  equation  (v  is  constant  and 
K  is  a  projection),  the  spectrum  in  the  spatially  uniform 
problem  (2.2)  is  the  single  Infinitely  degenerate  point  -v; 
in  the  initial  value  problem  corresponding  to  (2.6),  the 
spectrum  in  Pig.  3  becomes  a  vertical  line  through  -v;  for 
a  steady  flow  (2.7),  Pig. 4a  is  unchanged;  and  with  time 
separated  (2.8),  the  spectrum  in  Pig.  5  consists  of  two 
straight  lines  through  the  origin  (but  not  the  enclosed 
region).   The  spectral  resolution  (both  discrete  and 
continuum)  is  explicit  in  many  problems  for  the  relaxation 
model.   Several  steady  state  problems  have  been  solved 
with  this  technique  by  Cerclgnani  [ l4 ] ,  [153> 

In  the  cases  where  the  spectrum  falls  on  the  real 
axis,  the  identification  of  the  continuum  follows  from 
known  theorems,  and  where  it  is  complex  from  plausible  con- 
jectures (the  operators  are  not  self -adjoint ) .   The 
Indicated  perturbation  of  the  point  eigenvalues  is  a 


14.   C.  Cerclgnani,  Ann.  d.  Phys .  20,  219  (19^2). 

15-   C.  Cerclgnani  and  P.  Sernagiotto,  "Raylelgh's  Problem 
at  Low  Mach  Numbers  According  to  Kinetic  Theory,"  to 
appear  in  Proceedings  of  Pourth  International  Sym- 
posium on  Rarefied  Gas  Dynamics,  Toronto,  1964. 
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plausible  conject-ure  except  In  the  case  of  hard  spheres 
where  it  has  been  proved  for  sufficiently  small  k  [l6]. 
The  mathematical  proof  of  the  asymptotic  approximation 
of  Boltzmann  solutions  to  fluid  dynamics  is  Independent 
of  the  existence  of  such  perturbed  normal  modes  [9]. 

The  behavior  after  a  long  time  or  far  from  a  boundary 
Is  obtained  from  the  least  damped  part  of  the  spectrum. 
The  part  of  the  continuum  near  the  origin  always  comes 
from  large  values  of  4^.   In  Fig.  3a  the  dominant  effect 
is  clearly  given  by  the  point  eigenvalues.   In  Fig.  4 
the  ultimate  behavior  (linear  in  x)  is  given  by  the 
multiple  eigenvalue,  but  the  approach  to  this  asymptotic 
state  is  not  simple.   In  Fig.  5a  the  continuum  dominates 
very  far  from  the  wall,  but  this  contribution  is  small 
In  amplitude  (since  it  is  excited  by  the  tail  of  the 
distribution);  thus  the  more  highly  damped  eigenvalues 
dominate  for  a  while  if  co  is  not  too  large. 


16.   J.  A.  McLennan,  Physics  of  Fluids  8,  1580  (1965). 
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3.   Approach  to  the  Boundary 

There  is  a  nonuniform  approach  of  the  distribution 
function  to  its  value  at  the  wall  because  of  the  slow 
particles  which  have  a  very  small  mean  free  path.  We 
assume  here  that  a  solution  exists,  and  we  proceed  to 
find  its  properties.  Consider  a  steady  flow,  and  for 
simplicity  a  one -dimensional  problem.  The  steady  flow 
equation  (2.7)  can  be  written  in  the  integral  form 


f_^  =  f_^  exp(-vx/^^)  +  ^ 


X 

f 


exp(-v(x-y)/^^)K[f(y)]dy 

(3-1) 

for  particles  emerging  from  the  wall  (^-,  >  O),  and  a 
similar  expression  for  f_  (^-,  <  0).   Noting  the  identity 


X 

f 

V 


^J 


exp(-v(x-y)/^-j^)dy  =  1  -  exp(-vx/^^)  ,         (3-2) 


we  can  interpret  f  at  x  as  a  mean  value  of  f  and  the 
values  of  k/v  distributed  along  0  <  y  <  x  with  total 
weight  one.   The  relative  weights  assigned  to  f  and  K/v 
depend  on  x  and  |.   As  x  — >0  the  weight  concentrates  more 
on  f  ,  and  f  approaches  its  boundary  value  F  .   But  the 
approach  is  not  uniform  in  ^-,  ;  if  4-,  — >  0  as  x  — >  0  there 
remains  a  finite  deviation  of  f   from  f  .   This  nonuni- 
formlty  will  not  affect  any  quantity  integrated  over  | 
since  the  deviation  is  finite  and  is  concentrated  on  a 
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set  of  small  measure  as  x  — >  0.   Therefore,  if  we  assume 
that  f_  approaches  a  limiting  value  (whether  uniformly 
or  not),  the  integrated  expression  K[f]  will  approach  a 
definite  limiting  function  as  x  — >  0.   And  even  if 
f  =  (f+»f_)  is  discontinuous  at  ^-^   =  0,   K[f ]  will  be 
continuous.   For  purposes  of  the  limit  x  -^   0,  ^-,  -^   0 
we  need  only  consider 


k(e2.^)  =   11m  K[f]/v  (5.3) 

'1 


X — >o 
►o 


The  deviation  of  f  from  f   for  small  x  and  ^,  is  given  by 


f ^  -  f^  =  (k  -f^)[l  -  exp(-vx/43_)]  OA) 

where  v"  =  v{0  ,^^,i-^)  .      The  value  of  f ,  -  f,  at  any  (^^,i    ) 
depends  on  the  ratio  x/4-,  .   It  is  an  elementary  exercise 
to  verify  that 
+a 


[1  -exp(-v/x^-^)]d^^  =  2v  xlog  X  +  0(vx)      (3-5) 


J 


-a 


■X- 

In  a  plane  geometry  f  is  generally  discontinuous  at 
1^  =  0,  and  in  an  exterior  domain  across  a  line  of  sight 
cone.   But  in  a  convex  interior  domain  f  is  continuous  even 
in  free  flow  (e.g.,  with  a  varying  wall  temperature). 
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Thus  any  macroscopic  moment  of  f  -f   Is  0  {x   log  x), 
whereas  the  contribution  to  the  Integral  from  finite  ?,-, 

Is  only  0(x)- 

There  are  several  conclusions.  The  convergence 
of  f  to  Its  boundary  value  Is  nonuniform  In  a  small 
velocity  region  i-.         0  of  order  x.  Fig.  6.   There  Is  no 


Fig.  6 

reason  to  expect  the  approach  of  f_  to  be  nonuniform. 
The  approach  of  a  macroscopic  moment  to  Its  boundary 
value  Is  uniform,  x  log  x  — >  0,  but  the  derivative  Is 
unbounded.   This  Infinity  Is  In  addition  to  a  possible 
discontinuity  such  as  temperature  or  velocity  slip.  Fig.  ?. 
If  two  plates  at  different  temperature  or  velocity  are 
separated  by  a  distance  L  which  Is  small  compared  to  the 
mean  free  path,  the  Incoming  distribution  f  will  have  a 
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V 


Pig.  7 

disturbed  region  of  finite  width  0(L),  and  f  will  have 
a  bo-undary  layer  of  width  0(x)  as  x  — >  0.   A  macroscopic 
moment  will  have  Infinite  slope  at  each  wall  and  a  slope 
0(log  L)  at  the  mldplane . 

With  a  time  factor  1(d,  the  analysis  of  f   is  essen- 
tially the  same  except  that  v  is  replaced  by  v+loo.   The 
oscillatory  factor  exp( -ioox/l-,  )  yields  an  additional 
nonuniform  limit  with  very  rapid  oscillations  concentrated 
at  ^i  =  0.   A  typical  f  would  look  somewhat  as  shown  in 
Fig.  8.   If  00  »  1,  the  oscillatory  boundary  layer  with 
a  width   5^-,"^  cdx  dominates.   If  oj  «  1,  the  previous  boun- 
dary layer  of  thickness  bt-.^  x   dominates,  but  there  is  an 
oscillatory  sublayer  of  still  smaller  thickness  cux.   This 
complex  nonuniform  limiting  behavior  is  not  necessarily 
physically  Important,  but  it  must  be  taken  under  careful 
consideration  in  evaluating  mathematical  approximations. 


-  17  - 


Pig.  8 


4 .   Existence  In  the  Small 

Consider  a  one -dimensional  problem  between  two 
parallel  plates  for  the  separated  time  equation  (2.8). 
We  suppose  that  the  emerging  distribution  f  =  f   Is  given 
at  the  left  and  f_  =  T_   Is  given  at  the  right  (for  conver- 
sion of  these  results  to  a  more  conventional  boundary 
condition  which  relates  f  to  f_  at  a  wall,  see  Sec.  6). 
This  problem  is  essentially  the  same  one  that  has  been 
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treated  by  Willis  [l?];  the  presence  of  the  additional 
term  Icof  offers  no  increase  in  difficulty,  but  we  shall 
present  a  different  proof.   First  we  rewrite  the  Boltzmann 

equation  as  a  pair  of  Integral  equations 

X 


^+ 

=   f^exp 

(v+ia))x 

•i 

I      «i  - 

o 

f_ 

=   f"  exp 

(v+ico)  (x-L)' 

^^ 

L                  ^1 

J 

exp 


(v+io))  (x-y) 


Kdy 


X 


exp 


(v+lcp)  (x-y) 
^1    ■ 


(4.1) 
Kdy 


We  have  ?,.,  >  0  in  the  first  equation  and  ^,  <  0  in  the 
second.   For  a  steady  flow  we  merely  set  co  =  0.   We  shall 
wish  to  iterate  these  equations  for  small  L.   To  make 
estimates,  we  take  absolute  values 

X 


"  + 


< 


f + I exp 


vx 


^^ 


exp 


■_  v(x-y)' 
■  ^7^ 


K|dy 


(4.2) 


< 


exp 


v(L-xV 

~T7r 


+ 


exp 


v(y-x)' 

tItt 


K|dy 


X 


and  we  see  that  cd  disappears. 

We  cannot  drop  the  exponentials  to  simplify  the 


17.   D.  Roger  Willis,  "The  Effect  of  the  Molecular  Model 
on  Solutions  to  Linearized  Couette  Flow  with  Large 
Knudsen  Number,"  in  Rarefied  Gas  Dynamics,  ed .  L. 
Talbot,  Academic  Press  (19^1) • 
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estimates,  for  this  would  leave  the  non-lntegrable  factor 
1/4  .   We  know  that  K  Is  bounded  In  an  absolute  norm, 

|K[f]|  <  K^maxlfl  .  (4.3) 

But  we  cannot  use  this  simple  estimate  to  prove  convergence 
together  with  the  interpretation  of  (4.1)  as  a  mean  value 
of  f  and  K/v  because  k   is  larger  than  one.   We  cannot 
use  this  estimate  even  for  arbitrarily  small  L  because 
the  weight  [l  -  exp( -vx/^-,  )  ]  which  multiplies  K/v   is 
arbitrarily  close  to  one  if  ^,  is  small.   But  if  L  is 


smal 


1,  this  weight  is  small  for  most  values  of  ^,  .   We 


therefore  combine  a  maximum  and  Lp  norm  to  show  that  itera- 
tion of  (4.2)  about  the  second  term  on  the  right  converges 
for  small  L. 

We  use  two  properties  of  the  operator  K[f]  from  [7] 

(l+4^)^/^|K[f]|  <  K,max(l+^^)lf I  (4.4) 

|K[f]|  <  K2  [Jf^d^]^/^  (4.5) 

Both  K-,  and  Kp  are  larger  than  one.   Prom  the  first  rela- 
tion we  know  that  application  of  K  reduces  f  for  large  ^, 

2  "1/2 
specifically  for  (l+^  )  '   >  k.,  ;  but  K  can  increase  f  at 

smaller  values  of  ^ .   It  is  an  elementary  consequence  of 

(4.4)  and  (4.5)  that  if 
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fl  <  (1+^^)"^ 


1 


(4.6) 


A^  <  (Sk^k^)"^  =  a^ 


are  both  satisfied,  then 

|K[f]|  <  i  (1+^^)"^  .  (4.7) 

We  can,  in  addition,  find  a  value  L  such  that  for  |Kf| 
satisfying  (4.7), 

|[1  -exp(-vL/^^)]^  |K[f]|2d^  <  ^  a^  .      (4.8)* 
If  max  |f|  satisfies  (4.6),  we  say  it  has  norm  [l,a].   From 

X 

the  estimates  (4.7)  and  (4.8),  the  integral  terms  in  (4.2) 

have  norm  [1/2,  l/2a]  if  f  has  norm  [l,a].   The  iteration 

is  therefore  contracting  and  converges  uniformly  in  both 

the  maximum  and  Lp  norms. 

The  significance  of  the  values  of  the  various  constants 

depends  on  the  fact  that  the  Boltzmann  equation  has  been 

made  dlmenslonless,  distance  with  respect  to  mean  free  path, 

time  with  respect  to  mean  collision  time  [approximately,  by 

v(0)  =  1],  and  velocity  with  respect  to  mean  thermal  velocity. 

In  these  variables  the  constants  k  ,  k^ ,  k„  as  well  as  the 

o   1   2 

critical  value  of  L  are  all  of  order  one.   The  behavior  of 


* 

The  proof  of  this  inequality  is  not  difficult;  a  similar 
result  will  be  proved  as  the  Lemma  of  Sec.  5. 
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the  solution  In  co  Is,  of  course,  completely  lost  by  this 
method  of  estimation. 


5.   Existence  in  a  Half  Space 
We  consider  the  equation 

Icof  +  ^1  H  +  vf  =  K[f]  (5.1) 

in  the  half  space  x  >  0.   To  obtain  existence  for  unlimited 
X  we  cannot  afford  to  take  absolute  values  and  ignore  the 
oscillations  induced  by  00  (in  particular,  there  is  no  exis- 
tence theorem  for  steady  solutions,  00  =  0,  in  the  half 
space  X  >  0).   We  can  use  00  effectively  only  by  taking  a 
Fourier  or  Laplace  transform  in  x.   The  parameter  03  will 
reduce  the  amplitude  of  the  transformed  f  in  frequency 
space  except  at  certain  resonant  points.   We  cannot  afford 
to  use  a  maximum  norm  alone,  for  K  has  norm  larger  than 
one  and  cd  does  not  reduce  the  amplitude  uniformly.   An  Lp 
norm  will  be  reduced  by  cu  and,  as  in  the  previous  section, 
can  be  used  to  reduce  the  maximum  norm  as  well. 

Assuming  that  f  is  bounded  as  x  — >  <»,  we  can  rewrite 
(5.1)  as 
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f+  =  f^exp(-ax)+^  j  exp[-a(x-y)]K[f(e,y)]dy   (5-2) 

o 


f  = 


^ 


exp[a(y-x)]K[f(^,y)]dy 


(5.3) 


where 


a   =  J—   (v+Icd)  ; 
^1 


(5.M 


i-^   and  Re  a  are  positive  in  (5-2)  and  negative  In  (5-3) 
The  Incident  distribution  at  x  =  0  Is 


r 


f   =  - 


^ 


exp[a(y-x)]  K dy 


(5.5) 


o 


In  the  present  formulation  we  shall  assume  that  f   Is  given; 
f_  will  be  determined  by  the  condition  that  f  be  bounded. 
In  the  next  section  we  extend  this  formulation  to  Include 
more  conventional  boundary  conditions  relating  f  to  T  . 
We  Introduce  the  Laplace  transform 


f(^,s)  = 


e"^""  f(|,x)dx 


(5.6) 


i 


and  note  that  the  transform  of  K[f]  is  simply  K[f ] . 
Dropping  the  tilde,  we  have  the  transformed  equation 


f  = 


^-j^f  +  K[f  ] 
v+lco  +  i-.s 


(5.7) 
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a , 


We  would  obtain  the  same  result  by  a  Fourier  transform 
of  a  function  f  (x)  where  f  =  f  for  x  >  0  and  f  =  0 
for  X  <  0.   We  shall  use  s  for  the  complex  variable  and 
write  s  =  ik  with  real  k  when  s  is  to  be  restricted  to 
the  Imaginary  axis.   To  invert  the  transform  we  require 
that  f  be  analytic  for  Re(s)  >  0  and  f  — >  0  for  large  s. 
The  denominator  in  (5-7)  vanishes  If  -s  =  (v+ico)/^  - 
The  pole  is  in  the  left  half  plane  if  ^-,  >  0  and  in  the 
right  half  plane  if  ^^  <  0.   if  ^^  =  ^  =  0.  v  =  v(^^), 
this  pole  lies  exactly  on  the  curve  (Pig.  5)  which  was 

described  as  the  boundary  of  the  continuous  spectrum;  if 

2    2 

^2  +  ^^  >  0,  the  pole  is  to  the  left  of  the  spectral 

boundary  in  the  left  half  plane  and  to  the  right  of  this 
curve  on  the  right.   To  solve  the  problem  as  posed,  the 
singularities  in  the  right  half  plane  must  be  removed. 
To  this  end  we  set 


f  =  -  ^  K[f(|,-a)]  (5.8) 

^1 

for  ^.    <   0  and   split    (5-7)    into 

4.f,    +  K[f] 

f,    =  -^^ (5.9) 

v+icjo  +  s^i 

f      =   K[f(e,s)]   -K[f(^,-a)]  {b.\0) 

^-^(s+a) 

It  is  in  this  formulation,  (5.9)  and  (5. 10),  that  we  shall 
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solve  the  problem;  f  Is  assumed  to  be  known,  and  f  can 
be  computed  from  (5-8)  once  f  is  known. 

We  introduce  the  function  class  C  as  follows.   The 
domain  D_^  consists  of  all  real  ^^,    i    ,    ^.  >  0  together  with 
Re  s  >  0.   The  function  f_^(^,s)  is  in  C   if  (l+^^)f   is  con- 
tinuous and  bounded  in  D  ;  and  if  f   is  analytic  in  Re  s  >  0 
for  each  ^  (it  is  continuous  to  s  =  ik);  and  if  |f|  <  C/ | s | 
as  s  — >  00  for  each  ^  (C  need  not  be  uniform).   The  class  C 
is  similarly  defined  for  f_  and  D_  with  ^n  <  0.   We  say 
that  f  is  in  C  if  r_^   are  in  C_^.   There  is  no  implication 
that  f  is  continuous;  the  f  may  have  distinct  limits  as 
^1  ->to. 


f  =  f 


To  start  an  iteration  of  (5-9)  and  (5.10)  we  set 
(o) 


f|°^  =  ^1^+/  (v+icD+s^-^) 
(o)  _ 


(5.11) 


fy'  =  0 

We  remark  that  for  s  =  ik 

|v  +  i(cD  +  k^^)|  >  [1  +  (a3+k^^)^]^/2  >  1      (5.12) 

Taking  f  such  that  (l+^^)^  f  <  oo,  we  verify  that  f^°)  is 
in  C_|_.   It  is  continuous  in  (^,s),  it  is  bounded  on  the 
imaginary  axis,  and  it  is  analytic  and  decays  as  1/ | s |  in 
the  right  half  plane. 
Now  setting 
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^+        ~  ^+     +  V  +100  + s^-^  i^-j--^; 

V  +  loo  +  S^-j 

we  verify  that  all  Iterates  remain  In  C.   For  Induction  we 

fnl 
assume  this  property  for  f^  ' .      Prom  (^•'^)  we  see  that 

(1+4  )K^^'^  Is  bounded.   It  Is  continuous  and  analytic  In  s 

because  the  kernel  of  K  Is  continuous  and  bounded  (except 

for  an  Integrable  singularity  at  the  origin  [7])-   Using 

(5.12)  we  know  that  {l+i^)S^^^'    Is  bounded  on  s  =  Ik  and 

^(n+1)  __^  Q  ^3  |k|_^  00.   Therefore  fj^+^^  is  In  C_^.   We  may 

use  the  same  estimates  for  K^^'^(4,s)  and  K^"'^(4,-a)  In 
(5.14).   In  particular,  f^        '    Is  analytic  for  each  4;  the 
pole  Is  canceled  If  4-,  ^   0,    and  there  Is  no  pole  If  4,  =  0. 
Also  f ^"   '  Is  bounded  on  s  =  Ik  and  It  approaches  zero  as 


k 


oo;  for  4/0  this  Is  explicit,  and  for  |   =  0, 


K^"^(4,-a)  =  0  and  K^"^(4,lk)  ->  0  with  k.   It  remains  only 
to  establish  continuity  in  (s,4)  for  s  /^  oo  and  4-,'^  0. 
This  follows  from  uniform  convergence  on  s  =  Ik  as  4-1  — >  0. 
In  particular,  we  observe  no  difficulty  at  the  origin, 
s  =  0,  even  though  this  belongs  to  the  spectrum,  because 
this  part  of  the  spectrum  arises  for  large  4  where  all 
functions  are  uniformly  small. 

To  show  convergence  of  the  iterates  we  introduce  the 
same  norm  as  in  the  last  section,  viz.,  that  f  has  norm 
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[l,a]  If 

If|  <  (1+^^)"^ 

(5.15) 
[A^  <  a^ 

and  If  each  of  f  has  norm  [1,  a//2"].   We  set 

a  =  cVk^i1^2  (5-16) 

where  c  <  1  Is  a  positive  constant  to  be  fixed  later  [k, 
and  Kp  are  given  In  (^.•4-)  and  (4.5)],   We  now  prove  the 
Lemma:   Assume  that  f  has  norm  [l,a]  and  set 


g  =  [v  +  i(a)+k^-^)]  ^K[f]  .  (5.17) 


Then,  for  co  >  00   (depending  on  k  ,  Kp,  c),  g  has  norm  less 
than  [c  ,ca] . 

First  we  show  that  (1+^  )|g|  <  c.   Using  the  two 
bounds  (4.4)  and  (4.5)  In  turn,  we  have 

(1+^2) lgl  <  (i+e2)-V2  ^^ 

(l+^2)lg|  <  {l+t^)K^a 

2     2   2 
Setting  1+4   =  K-j/c  ,  the  first  estimate  Is  less  than  c 

for  ^  y   i     while  the  second  Is  less  than  c  for  ^  <  ^  .   To 

f     2  2  2 

complete  the  lemma  we  must  show  that  J  g  d^  <  c  a  .   We 

have,  using  (4.4)  and  (5-12) 
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00  00 

=  ttkM  [1+  (oo+k^-L)^]"^d^-j_  J  [1  +  u  +^^]"^du 


-00  O 


This  can  be  Integrated  explicitly,  but  It  Is  easier  to 

estimate.   For  the  limited  range  of  Integration  where 

2     2 
|co  +  k:^-,|  >  \x,   writing  [  1  +  (cu  +  k^.,  )  ]  >  |J-  we  obtain  the 

upper  bound 

1   2  2/2 

j^  IT    K^/^i   . 

For  the  remainder,  |Go  +  k^-,|  <  i^,  we  write  [  1  +  (a)+ k^-,  )  ]>1 

2  2   2 

and  {1  +i-.)   >  [1  +  (co-l-t)  /k  ]  to  obtain  (after  maximizing 

with  respect  to  k) 

Sir  ^2    H 

^  ^1  Ico  -m      • 

Together, 

2  r  TT^   1  ,  Stt    ^l 


g  d^  <  K^  L  T^  -^  +  25  |a)  -^l|  -I 


Setting 

m     =  II  +  32ti^/25Tr 


2    2  2/„  2  2 
H-  =  TT  K-,/2a  c   , 


(5.18) 


/■   2       2  2 

we  have  I  g  d^  <  a  c  for  cu  >  cu  ,  proving  the  lemma, 
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We  are  now  ready  to  show  convergence  of  the  Iteration 
for  o)  >  00^.   If  fj"^  -  f[^~^)  has  norm  [1,  a//2  ]  then 
fjn+1)  _  ^|n)  ^^^  ^^^^  ^^^    ca//2  ]  .   Similarly,  the  norm 
[1,  a//2  ]  for  f}^^    -   f(""^)  Induces  a  norm  [2c,  /2"  ca] 
for  f _    '^  -  f ^  '  .   The  iteration  is  contracting  if  we 
take  c  <  1/2  and  take  f_J^  '    of  finite  norm.   It  is,  of 
course,  sufficient  to  make  these  estimates  on  the  imaginary 
axis  s  =  Ik  since  all  the  Iterates  are  analytic. 

By  the  uniform  convergence  we  have  proved  the  exis- 
tence of  a  solution  f(^,s)  which  is  in  C   The  Inverse 
transform  exists  and  yields  a  bounded  solution  in  x.   The 
iterations  also  converge  uniformly  to  the  incident  f_.   And 
since  the  convergence  factor  increases  with  o),  the  function 
f_  goes  to  zero  for  large  oo.   In  other  words,  no  particles 
are  scattered  back  when  oo  is  large.   Also  the  first  term 


f_;  '    is  dominant  in  some  sense  (which  we  shall  examine  in 


Sec.  7)  for  large  co.   The  iteration  of  the  Laplace  transform 
f(^,s)  corresponds  exactly  to  Iteration  of  the  Integral 
equations  (5-2)  and  (5.3)  for  f(^,x). 

A  better  bound  on  f   such  as  f  ~  (l  +  ^  )    will  yield 
corresponding  estimates  for  the  solution.   The  decay  as 
X  — >  00  can  be  estimated  by  integrating  the  inverse  Laplace 
transform  by  parts.   If  f  decays  algebraically  for  large  ^, 
then  f  will  decay  algebraically  in  x  [f(^,s)  has  only  a 
finite  number  of  derivatives  at  s  =  0].   If  f  has  all 
moments  (e.g.,  Maxwellian) ,  then  f  will  decay  transcendentally 
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in  X. 

An  existence  theorem  for  the  nonlinear  Boltzmann 
equation  can  be  carried  out  in  almost  the  same  way,  using 
the  collision  term  estimates  of  Ref,  [8].   The  result  is 
surprisingly  uninteresting.   In  ordinary  gas  dynamics  a 
linearization  becomes  Increasingly  worse  at  high  frequencies 
(the  common  linearization  takes  derivatives  to  be  uniformly 
small).    But  the  correct  kinetic  approach  shows  that  at 
very  high  frequencies,  where  the  fluid  analysis  is  no 
longer  valid,  the  nonlinearity  becomes  progressively  less 
important.   The  reason  is  that  in  the  Boltzmann  equation 
the  nonlinearity  resides  only  in  the  collision  term, 
whereas  it  is  the  strictly  linear  streaming  term  which 
dominates  at  high  frequency.   In  other  words,  the  solution 
decays  to  linearity  in  much  less  than  a  mean  free  path 
before  collisions  play  a  significant  role.   This  is  an 
important  point  with  regard  to  experiment,  because  high 
frequency  experiments  are  inevitably  operated  at  a  high 
excitation  level  in  order  to  be  able  to  detect  the  strongly 
damped  waves. 
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6.   Boundary  Conditions 

The  problem  that  we  have  solved  Is  artificial  in 

that  the  emitted  distribution  f,  was  assumed  to  be  given 
outright.   The  more  conventional  formulation  is  to  give  f" 
In  terms  of  f_.   The  general  form  of  such  a  boundary  con- 
dition is 

f^  =  f Q  +  B[f_]  (6.1) 

where  f   is  given  and  B  is  a  bounded  linear  operator;  we 
shall  give  B  explicitly  for  specular  and  diffuse  reflection 
below.   But  the  solution  obtained  in  the  last  section  gives 
f_  as  a  bounded  linear  operator  on  f  , 

f_  =  R[f^]  .  (6.2) 

This  operator  is  obtained  explicitly  In  the  iteration  as  a 
geometrically  convergent  sequence,  and  the  norm  of  R  goes 
to  zero  for  large  co.   The  linear  functional  equation 

f^  =  f Q  +  BR[  f^]  (6.3) 

has  a  convergent  Neumann  series  for  sufficiently  large  cd, 
viz.,  for  which  the  norm  of  BR  is  smaller  than  one.   There- 
fore f  can  be  found  and  the  existence  theorem  applied  to 
the  general  boundary  condition  (6.1). 

In  a  specific  problem  it  would  be  more  expedient  to 
iterate  (6.1)  and  the  Boltzmann  equation  together,  taking 
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T  =  f  as  the  Initial  boundary  condition  and  improving 
f  and  f  in  turn.   The  dominant  term  for  very  large  co 
would  be  the  combined  first  Iterate, 


f^  =  fQexp[-(v  +  loo)x/e^] 

(6.4) 
f_  =  0 

Both  the  collision  operator  K  and  the  boundary  operator  B 
become  unessential  in  the  limit  of  large  oo. 

With  a  diffuse  reflection  boundary  condition  the  full 
nonlinear  boundary  distribution  is 


+    (27rRT)^/'^ 


Ilzlil 


2 


2RT 


(6.5) 


where  u(t)  represents  the  oscillating  wall  (both  normal 
and  tangential),  T(t)  can  also  oscillate,  and  p(t)  must  be 
left  open  to  satisfy  conservation  of  mass.   We  linearize 
according  to  the  scheme 

/  -,    ,    icot  \ 
P  =  Po(l  +Pie    ) 

T  =  T  (1+T,e^^^)  (6.6) 

o  ^     1     ' 

u  =  (RT  )  '^   u.e 
^   o'      1 

_     -1/2 
and  introduce  the  dimensionless  velocity  ^  =  ^  (RT  )  '^ 

and  the  normalization 
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P  =  F  (1  +0)^/^  f ) 


^  = -^j^   exp(-  ^  ^  ) 


(6.7) 


(2.)"^ 


which  are  implicit  in  the  derivation  of  the  dimenslonless 
linear  Boltzmann  equation  which  we  have  been  using  (for 
details  see  [10]).   We  have 

-   l/2r,  ^^  ,1  ,2   3, 


f 


+ 


=  CD  '  [p-L  +T^(f  r  -f)  +  u^^^  +  M^^^-\    ,  (6.8) 


where  u-,  is  the  normal  component  and  Up  the  tangential 
component  (taken  for  simplicity  in  the  direction  ^p). 
The  conservation  of  mass, 

(6.9) 


I   ^jf^oD  d^  +    ^^_^   d|  =  0 


^^>0  ^^0 


serves  to  eliminate  p-j^  and  gives  the  boundary  condition 


TT       1/2 

f   =  CO  '^ 


-  (277X1))^/^ 


(6.10) 


e^f_a)  d^ 


The  first  term  [f^  in  (6.1)]  is  linear  in  the  given  oscil- 
lation amplitudes  u-j^,  u^ ,  and  T-^   and  is  a  polynomial  in  i 

1/2 
times  03^  .   The  operator  B  in  (6.1)  is  in  this  case  a 

simple  projection  (first  moment  of  f  ) .   For  a  simple  normal 
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oscillation  of  the  wall  we  set  u^  =  T-,  =  0,  for  a  simple 
shear  u^  =  T,  =  0,  and  for  a  stationary  but  heated  wall 
U-,  =  u„  =  0.   Since  the  formulation  Is  linear,  each 
elementary  problem  can  be  solved  separately. 

A  similar  linearization  of  the  specular  reflection 
boundary  condition  of  a  normally  oscillating  wall  (neither 
shear  nor  temperature  variation  can  couple  to  the  gas  with 
specular  reflection)  yields  the  boundary  condition 

I^ii^)   =   cD^/2|^u^  +  f_(-^i)  (6.11) 

The  operator  B  in  this  case  is  the  identity. 
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7-   High  Frequency  Behavior 

The  dominant  feature  at  high  od  Is  that  f  Is  rapidly 
oscillating  In  4  as  well  as  In  x.  Pig.  9.   The  first 


Pig.  9 


iterate  is  given  by 


f_^  -  fQexp[-(v  +  ia3)x/e^] 


f  =  0 


(7.1) 


where,  as  In  {6A),   we  have  taken  the  leading  term  in  the 
boundary  condition  as  well  as  in  the  Boltzmann  equation. 
If  we  increase  cd  at  fixed  x  and  ^,  f  oscillates  but  with 
fixed  amplitude  f ^exp( -vx/^^ ) .   We  shall  verify  that  at 
any  fixed  x  and  4,  i^  ^   0,  the  error  in  taking  the  leading 
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term  (7-1)  Is  small  compared  to  this  amplitude  factor 
f  exp( -vx/^i  ) .   Any  moment  or  other  Integrated  quantity- 
over  f   Is  small;  algebraically  as  co   If  f  decays 
algebraically  for  large  |,  and  transcendentally  small  in 
CD  if  f  decays  transcendentally  in  ^  (e.g.,  Maxwellian). 
Thus  K[f  ]  and  the  consequent  error  in  the  formula  (7.1) 
are  small  compared  to  the  term  which  is  kept, 

lim  |f_^  -  f^exp[ -(v+ ia3)x/^-j^]  I  -  0 

(JO— >o 

lim  If  I  =  |f^lexp(-vx/^^)  (7.2) 

CL)-*0 

lim  |f_|  -  0 

CD— >0 

This  result  is  interesting  in  that  x  can  be  several 
mean  free  paths.   Collisions  do  have  an  effect  in  cutting 
down  the  boundary  emission  f  by  the  factor  exp( -vx/^-,  ) , 
but  there  is  no  compensating  term  from  the  creation 
operator  K[f]  to  this  order.   The  reason  is  simply  that  vf 
is  singular  and  K[f]  is  regular;  any  averaging  over  f  reduces 
the  amplitude  because  of  the  rapid  oscillations. 

The  answer  is  much  more  subtle  if,  instead  of  the 
error  in  f ,  we  look  at  the  error  in  a  moment  of  f .   In 
particular,  it  is  a  moment  which  is  usually  the  experimental 
observable.   For  a  Maxwellian  f  ,  any  moment  will  be  trans- 
cendentally small  in  co.   It  is  plausible  that  the  next 
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Iteration,  which  first  Involves  evaluation  of  K[f]  (compar- 
able to  evaluating  a  moment)  and  then  another  average  of  K 
with  a  highly  oscillating  Integrand  [cf.  (4.1)],  will  give 
a  contribution  smaller  than  the  leading  term.   For  an 
algebraic  decay  of  f^  this  can  be  demonstrated.   But  for 
the  much  more  relevant  exponentially  decaying  f  ,  this 
estimate  presents  great  difficulty.   The  reason  Is  that 
the  leading  term  (7-1)  cannot  approximate  the  actual  solution 
uniformly  for  small  ^-^^   I/odx.   Even  for  arbitrarily  small  x 
we  have  seen  that  f_^  does  not  uniformly  approximate  f 
(Sec.  3).   If  the  moment  is  transcendentally  small  compared 
to  the  fluctuating  amplitude  of  f,  we  cannot  eliminate  the 
small  region  ^-^ —  l/oox  from  the  integration  without  estimating 
it  accurately.   Although  this  is  a  gap  in  the  adoption  of  the 
leading  term  (7-1)  for  integration  purposes,  the  area  of 
doubt  in  the  estimate,  viz.,  small  ^,,  is  where  f  (?,) 
oscillates  most  rapidly.  Fig.  9;  it  is  quite  plausible  that 
this  difficulty  is  technical  rather  than  actual. 

Summarizing,  the  amplitude  of  f^°^  remains  finite  as 
CO  — >0  and  is  clearly  dominant  asymptotically;  the  estimate 
is  more  subtle  for  a  moment  of  f  but  is  very  plausible. 

We  shall  now  accept  the  accuracy  of  (7.1)  and  evaluate 
its  moments  (e.g.,  density  and  pressure)  for  possible 
experimental  comparison.   For  either  diffuse  or  specular 
reflection,  the  function  f  has  the  form 
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_ 1  ?2  1  >2 

IT  ^  P    IT  ^ 

f  -  ^.  e  ^  or  r  e  ^  (7.3) 

o    1 


for  large  t,   depending  on  whether  the  wall  Is  moving  or 

is  given  a  temperature  oscillation.   Before  integrating, 

2 
f  is  multiplied  by  1  or  ^   depending  on  whether  density 

or  pressure  is  measured.   Thus  the  function  to  be  evaluated 

has  the  form 


Q 


=  I  ^''exp[-  J  i^   -vx/^3_  -icox/^-^]d^       (7.^) 


where  n  takes  values  between  1  and  4.   The  collision 
frequency  for  large  ^  has  the  behavior 


•  a 


v(l)-v^^"  (7.5) 

where  v-,  is  of  order  one  and  0<a<l  (a  =  0  for 
Maxwellian  fifth  power  and  a  =  1  for  hard  spheres). 

Any  point  of  steepest  descent  is  found  at  ^„   -   ^^  -   0; 
thus  we  can  consider  the  integration  (7-3)  in  one  dimension 
and  ignore  the  distinction  between  |  and  ^-,  .   A  point  of 
steepest  descent  is  a  stationary  point  of 

-,      o    v^x 
-1    fi    1  e2      1     ia)X  fry   r\ 

n  log  ^  -  If  ^   -  -^  -  -r  '  ^^-^^ 

The  dominant  contribution  for  large  co  is  free  flow,  given 
by  the  second  and  fourth  terms  (see  also  [6]), 

il  =   2icDX  .  (7-7) 


o 
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The  correct  critical  points  are 


i^  =  (2a3x)i/^(|  it^yj) 


(7.8) 


and  the  path  of  integration  for  i   Is  as  shown  in  Pig.  10. 


Pig.  10 

The  dominant  asymptotic  expression  for  Q  is  oscillatory 
in  oox  but  with  an  amplitude  given  by 


log  Q I  (2a)x)^/^ 


The  wavelength  of  the  oscillation  of  Q  is  locally 


(7.9) 


A 


^TT  /2x\ 


1/3 


(7.10) 


which  is  short,  but  slowly  lengthens  with  x. 

If  X  is  considered  to  be  fixed  as  oo  Increases,  the 
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next  order  term  taken  from  (7-5)  is  n  log  X  and  the 
effect  of  collisions  only  comes  In  at  third  order.   To 
see  the  effect  of  collisions  more  simply,  we  let  x  grow 
as  (J''      which  brings  in  v(^)  at  second  order.   We  obtain 

log  Q^-  2  (2oox)^/^  -  v,x(2a)x)"^^""^/^cos[5'n-(l-a)/6] 

O  -L 

(7.11) 
(valid  for  any  scaling  of  x  which  makes  vx/l  dominate 
log  ^)  and 

log  Q--|  (2)2/^co  -  V3_2-^^-")/^cos[57r(l-a)/6]cD°^2 

(7.12) 

1/2 
for  X  =  CD  '  .   An  experimental  confirmation  would  have  to 

detect  two  terms  in  a  plot  of  log  Q.   The  first,  linear 

in  CD,  would  be  a  confirmation  of  the  Maxwellian  distribu- 

a/2 
tion.   The  second,  proportional  to  ao  '   would  give  a 

measurement  of  the  force  law  exponent  a  and  coefficient  v-^ 

for  large  velocity  ^ . 

Taking  a  second  term  in  the  expansion  with  x  fixed, 

we  have 

log  Q I  (2cDx)^/^  +  I  n  log  CD  .        (7-13) 

Observation  of  the  second  term  with  this  scaling  would 
show  the  deviation  from  a  Maxwellian  evidenced  by  the 
number  n. 
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In  summary,  any  direct  measurement  of  the  distribution 
function  would  give  a  nontrlvlal  lowest  order  result  involving 
the  entire  function  v(^)  and  the  boundary  condition  through 
the  expected  variation  f^exp( -vx/^-^) .   Measurement  of  a 
moment  of  f  gives  information  about  the  boundary  condition 

if  we  keep  x  =  constant  and  information  about  v(^)  for  large  ^ 

1/2 
if  we  scale  x  /^  cu  ^  ,   All  measurements  will  be  amplitudes  of 

rapidly  oscillating  fluctuations  in  x,  i,   cd  as  well  as  t. 

We  do  not  attempt  to  compare  these  results  with  experi- 
ment.  The  lowest  order  free  flow  expression  (7-9)  has 
apparently  been  observed  experimentally  [l8,  19,  6].   Since 
the  dimensional  x  and  cd  vary  as  p  (pressure)  and  l/p 
respectively,  in  a  given  experiment  p  must  be  varied  and  in 
addition  either  x  or  cd  must  be  varied  in  order  to  be  able 
to  observe  the  second  term  in  (7. 12). 

We  remark  that  an  approximation  such  as, that  made  by 
Kahn  and  Mintzner  [6],  but  with  our  f^°^  as  the  starting 
point  rather  than  free  flow,  might  yield  very  accurate  high 
frequency  data. 


18..  E.  Meyer  and  G.  Sessler,  Zeitschrift  f.  Physik  1^9,  15 
(1957).  

19.   G.  Maidanik  and  M.  Heckl,  Physics  of  Fluids  8,  266  (I965) 
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